Dissipative dynamics of a kink state in a Bose-condensed gas. 
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We develop a theory of dissipative dynamics of a kink state in a finite-temperature Bose-condensed 
gas. We find that due to the interaction with the thermal cloud the kink state accelerates towards 
the velocity of sound and continuously transforms to the ground-state condensate. We calculate the 
life-time of a kink state in a trapped gas and discuss possible experimental implications. 



I. INTRODUCTION 

The recent successful experiments on Bose-Einstein 
condensation (BEC) in trapped clouds of alkali atoms 
have boosted an interest in the physics of ultra- 
cold gases 0] . One of the challenging goals is to create 
and study macroscopically (topologically) excited Bose- 
condensed states, such as 2D or 3D vortices and their 
ID analogies, the so-called kinks. These excited states 
attract a wide attention because they behave as particle- 
like objects (solitons) and thus bring in the analogies with 
high-energy physics, where elementary particles are sug- 
gested to be solitons of fundamental fields. 

The simplest example of an excited condensate state 
is a kink- wise state in a cylindrical harmonic trap || [?|. 
This state has a macroscopic wavefunction with one 
nodal plane perpendicular to the symmetry axis and rep- 
resents a non-linear wave of matter. The investigation of 
these waves is especially interesting for understanding to 
what extent they are similar to light waves in non-linear 
optics and what is the difference related, in particular, 
to the interparticle interaction in a matter wave. 

Fundamental limitations for creating and studying 
kink states concern their stability. First of all, a kink 
state is always thermodynamically unstable, since its en- 
ergy is larger than the energy of the ground-state Bose 
condensate. In contrast to the well-known kink solution 
in ID in three dimensions the kink state is also char- 
acterized by the presence of a transverse dynamical in- 
stability, i.e. the instability of small-amplitude normal 
modes of the motion parallel to the nodal plane. In order 
to suppress this instability in a trap, one should strongly 
confine the radial (transverse) motion and turn to an ef- 
fectively ID case by making the radial frequency larger 
than the mean- field interparticle interaction . It is this 
condition that makes a kink- wise condensate stable in the 
limit of zero temperature. 

The thermodynamic instability will lead to decay of the 



kink state only in the presence of dissipative processes. 
Those are related, for example, to the interaction of a 
kink with a thermal cloud. In this Letter we analyze the 
dissipative dynamics of a kink state in a Bose condensate 
at finite temperatures. We show that the thermal excita- 
tions scattered by a kink stimulate its diffusive motion, 
with an increasing mean velocity, and induce the dissipa- 
tive force proportional to the velocity of the kink. After 
a short time the dissipative force takes over the diffusion 
and the kink steadily accelerates towards the velocity of 
sound. In the course of accelerating the kink state con- 
tinuously transforms to the ground-state condensate. 

Importantly, this process is insensitive to the size of the 
condensate in the direction perpendicular to the nodal 
plane. The key features remain the same as in an in- 
finitely long condensate. This is fundamentally different 
from the dissipative dynamics of a vortex state ||, al- 
though in some sense the interaction with the thermal 
excitations is similar. The vortex state has a topological 
charge (circulation) and in an infinitely large system can 
not continuously transform to the ground-state conden- 
sate. The vortex can loose the charge and, hence, decay 
only when approaching the border of the condensate. 



II. MEAN-FIELD DYNAMICS OF A KINK STATE 

Wc first briefly outline the mean-field (frictionless) dy- 
namics of a kink state as a solution of the ID Gross- 
Pitaevskii equation (GPE) in a homogeneous condensate 
of density n Eo] : 



ih 



9£o 
dt 



2m dx 2 



(1) 
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where g = 4:irh 2 a/m, a is the s-wave scattering length, 
m the atom mass, and y, = gno is the chemical potential. 
The condensate wavefunction \&o exp (-ifj,t/h) in Eq.(§) 
is characterized by a single distance scale, the correlation 
length Iq = h/ yjmnog. The validity of Eq.(|l]) for describ- 
ing spatial variations of requires Iq much larger than 
the mean interparticle separation along the x axis. The 
corresponding criterion can be written as 



n a Slo > I, 



(2) 
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where S is the transverse cross section of the condensate. 
For the kink state with the y, z nodal plane moving with 
a constant velocity q along the x axis, Eq. (|l|) gives |TT|j[i[ 

Mx) = ^(i± + ^t^^£j^. (3) 

Here q is the position of the nodal plane, cs — \fn®g/m 
is the speed of sound, and it is assumed that \q\ < eg. 

The energy of the kink state (|^) can be represented in 
the form (see jfjj) 

H(q) = ~Mc| (l - J) ; M = A(n Sl )m. (4) 

For the velocity of the nodal plane, \q\ <C c$, we find 

H(q) = hlc 2 - X -Mq\ (5) 

On the basis of Eq.(j|), the quantity — M can be treated 
as a negative mass of the kink. The criterion (||) ensures 
that M 3> m, i.e. the kink as a particle-like object is 
much heavier than a gas particle. The kinetic energy 
term in Eq.(|5|) is negative, and the energy of the kink de- 
creases with increasing velocity q. This means that if the 
dissipation decreases the kink energy, then the kink ac- 
celerates, gradually reaching the velocity of sound. Once 
this happens, the wavefunction (||) becomes nothing else 
than the wavefunction of a motionless ground-state con- 
densate. Accordingly, the time scale on which the kink 
accelerates towards the velocity of sound can be regarded 
as a life-time of a kink- wise state. 

The equations of motion for the kink as a particle-like 
object follow from the Hamiltonian H (0) expressed in 
terms of momentum p. The relation between p and q is 
obtained straightforwardly from the Hamiltonian equa- 
tion p = dH/dq: For slow kinks we have p = —Mq 
and H(p) = Mc|/3 -p 2 /2M. Eq.(g) demonstrates an 
important kinematics property of a kink state. As the 
time derivative p — — M(l — q 2 / 'c 2 s ) l / 2 q, the quantity 
— M(l — q 2 /c^) 1 / 2 can be viewed as a kinematic mass 
of the kink. This kinematic mass decreases when q ap- 
proaches the sound velocity cs, and it becomes easier to 
accelerate the kink. This is opposite to the case of a 
relativistic particle, where the kinematic mass increases 
with q, and an infinite force is required to accelerate the 
particle beyond the velocity of light. Here, the kink can 
ultimately reach the velocity of sound and disappear. 

III. KINETIC EQUATION FOR DISSIPATIVE 
EVOLUTION 

We now turn to the analysis of dissipation in the pres- 
ence of a thermal cloud and treat the kink as a particle- 
like object. The motion of the kink in a dissipative envi- 
ronment is not deterministic and is characterized by the 



probability density F(t,p, q) of having the momentum p 
and coordinate q at time t. The distribution function F 
satisfies the kinetic equation (see, e.g. Jl2[) 

dF d fdH\ d fdH\ T „ 

lH~Tp{^ F ) + d- q {-dp- F )= I[F] - (6) 

The colhsional integral I[F] in the r.h.s. of Eq.(^jJ) de- 
scribes the interaction of the kink with thermal excita- 
tions and, hence, is responsible for the dissipation. In the 
absence of dissipation (I[F] = 0) the l.h.s describes the 
motion of the kink, governed by the Hamiltonian equa- 
tions p = —dH/dq, q — dH/dp. 

The colhsional integral accounts for the change of the 
kink momentum and energy due to scattering of thermal 
excitations from a moving kink. As the kink is much 
heavier than the gas particles (M ^> to), one can expect 
that in each scattering event the relative change of the 
kink momentum p is small. Then, similarly to the case 
of motion of a heavy particle in a gas of light particles, 
the colhsional integral takes the Fokkcr-Plank form [Q 

The first term in the l.h.s. of Eq.(Q) describes the mo- 
bility in the momentum space. The transport coefficient 
A represents the mean momentum transferred per unit 
time from the thermal excitations to the moving kink, 
and hence is equal to the friction force acting on the 
kink. The second term is responsible for diffusion in the 
momentum space, with B being the diffusion coefficient. 

We will assume that the characteristic energies of 
transverse excited modes greatly exceed the gas tempera- 
ture T . Accordingly, only longitudinal (and in this sense 
ID) excitations are present in the system and interact 
with the kink. In the reference frame, where the kink is 
at rest, the reflection of an excitation from the kink pre- 
dominantly changes the excitation wavevector k to —k. 
This automatically follows from the condition M ^> to. 
Then, as in each scattering event the momentum trans- 
ferred from the excitations to the kink is equal to 2hk, 
the equations for the transport coefficients read 

where e& = [E\ + 2yxE/ c ) 1 / 2 is the Bogolyubov energy of 
the excitations, Ek = Ti k 2 /2m, R(k) is the reflection co- 
efficient of an excitation with momentum hk, and N(ek) 
are equilibrium occupation numbers for the excitations. 
The quantity in curly brackets represents the scattering 
rate for the excitations with wavevectors in the interval 
from k to k + dk. For subsonic kinks (\q\ <ti cs), we have 
hkq <C £fc. Then, in the first integral the occupation 
numbers N{e^ — hkq) have to be expanded in powers of q 
up to the linear term, which gives the mobility coefficient 
A = —aq. In the second integral one can put q = 0, thus 
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omitting an inessential dependence of the diffusion coef- 
ficient B on q. As a result, for the (positive) transport 
coefficients a and B we can write 



-2h 



B = 2h 



° d A R{k) d J±e d JL 



(8) 
(9) 



The ratio B/a = e* , where e* ~ min{T, /i} is a char- 
acteristic energy of the excitations which give the main 
contribution to the transport coefficients (see below). 



IV. DIFFUSION- 



AND MOBILITY-DOMINATED 
REGIMES 



We first analyze qualitatively the kinetic equation (g) . 
For an axially homogeneous condensate dF/dq — 0, and 
using Eqs.(|[), (0) we obtain the equation describing the 
evolution of a subsonic kink (p = —Mq): 



dF 
~dt 



d 



Mdq 



aqF 



B dF 



(10) 



At small times t the diffusion dominates over the mobil- 
ity and leads to increasing mean square velocity: (g 2 )~ 
Bt/M. Accordingly, the ratio of the mobility to diffusion 
term in Eq.^) also increases. Once (q 2 ) reaches a crit- 
ical value q 2r ^B /Ma~e* / M , the mobility takes over the 
diffusion. This happens on a time scale 



t ~ td = Ma 



(11) 



Owing to the condition M ^> m, the velocity go *C cs, 
i.e. the kink remains subsonic when the mobility starts to 
dominate over the diffusion. The condition M 3> m also 
justifies the use of the Fokker-Plank approach: The kink 
momentum Mqo can be written as (M/i/me*) 1 ' 2 e*/cs 
and it greatly exceeds the momentum 2hk ~ e*/cs trans- 
ferred to the kink in a single scattering event. 

At much larger times, t > td, the diffusion spread- 
ing of the kink velocity is negligible compared to the 
mobility-induced increase of q. As the force A acting on 
the kink is proportional to the kink velocity q, the latter 
grows exponentially with time. Omitting the diffusion 
term, Eq.(|^) is reduced to 

f)F r) 

and coincides with the Liouville equation for particles 
moving in accordance with the law p = A. Then, with 
A = —aq and p = —Mq, wc find q = q(0) exp (£/t£>) or 
t = td log (q/q(0)). From Eq.(||) one can see that 



IF 



-aq 



(12) 



Eq. (|12() clearly shows that the decrease of the kink energy 
is equal to the energy losses due to the friction. 

The logarithmic dependence of t on q ensures that the 
dissipative dynamics of a kink state has a bottleneck at q 
significantly smaller than cs- With logarithmic accuracy, 
this result is also valid at q ~ cs, unless q is very close 
to cs in which case the kink state is practically indistin- 
guishable from the ground-state condensate. Therefore, 
a characteristic time r at which the kink acquires the 
sound velocity and disappears can be found from the ob- 
tained dependence t{q), with q — eg: 



T D \og(c s /q(0))- 



(13) 



The life-time r of the kink is much larger than td and is 
insensitive to a precise value of the initial velocity. This 
allows one to put q(0) ~ q Q ~ (e^/Af) 1 / 2 . 



V. REFLECTION OF EXCITATIONS FROM A 
KINK 

For finding the time td and velocity go (energy e*) one 
should calculate the transport coefficients a and B. This 
requires us to solve a ID scattering problem and find the 
reflection coefficient R(k) of an excitation from a static 
kink. Following the Bogolyubov approach, we represent 
the non-condensed part of the field operator of atoms 



, b k are the annihila- 



as * = J2k( u kh - v k b\), where h 

tion/creation operators of elementary excitations. The 
wavefunctions f£ = u^ivk of the excitations satisfy the 
Bogolyubov-de Gennes equations 

£fe/ * 7 = ~&& f * +5l*oW| 2 /+ - M/ fe + , (W) 



2m dx 2Jk 



^ft = — — h + 3.9l*o(x)| 2 / fc - - M / fe - (15) 



The condensate wavefunction \E"o is given by Eq. ([|) with 
g = 0, and l^ol 2 is characterized by the presence of a 
density dip of spatial size lo, associated with the kink. 

For the phonon branch of the excitation spectrum 
(en <C n) we have solved Eqs.(|l4|),(p"5|), relying on the 
presence of the small parameter The method is 

based on matching the wavefunctions of free motion of 
Bogolyubov excitations in the range of distances x from 
the kink, where Iq -C |x| -C fc _1 , with the excitation 
wavefunctions obtained from Eqs (|l4|),([l^) by the per- 
turbation expansion in powers of e^/fj,. 

The Bogolyubov-de Gennes equations (jlj) , ( |l5| ) have 4 
fundamental solutions corresponding to = 0, k = : 



ft = tanh^/Zo); /r = 0, 

!t = 0; fz = cosh~ 2 (x/l ), 

ft = (x/Zo) tanh(x/Z ) - 1; f^ 

ft = 0; / 4 - = sinh(2x/Z ). 



0, 



(16) 



The modes ff, are odd, and the modes f%, ff are 
even with respect to changing the sign of x. For finite 
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satisfying the condition e k <C \i {klo <C 1), at \x\ <C k^ 1 
each particular solution can be written in the form of ex- 
pansion in powers oie k //j, around one of the fundamental 
modes. The coefficients of the expansion are obtained di- 
rectly from Eqs.([r^),(|l5|). For finding the solution tend- 
ing to ft (or to f^) as eft — > 0, we substitute ft (or f£) 
into Eq.jjlq) and, performing the integration, obtain the 
contribution to fZ, proportional to e k /fx. The function 
fZ is then substituted into Eq.(|l4]), which allows us to 
find the contribution to ft, proportional to (e k / fi) 2 , and 
so on. A similar procedure is used for obtaining the solu- 
tions tending to f 2 and f± as e k ~* 0. For the solutions 
ft and , which in the limit of zero correspond to 
/ x and f 3 , it is sufficient to keep the terms independent 
of and linear in e k /fx . Then we find 



/ 3 + ; /s~ fc = {cosh 2 (x/; ) + 2/. 

+ [(xAo) 2 + C 3 ]/ 2 -}, 



(17) 
(18) 



where the fundamental modes ff are given by Eqs.([16|), 
and C3 is the integration constant. The fundamental 
mode f% decays exponentially at large \x\. Hence, for 
the particular solution going into /^as e k — > we should 
also retain quadratic terms in the expansion. This gives 

/ 2 + fe = (WM)[(C lfc -2)/ 3 +-l]; (19) 
f-2k = f-2 + ( e k/2fi) 2 {Cifccosh 2 (a;//o) 

+ (C lk - 2) (2/ 3 + + 1 + [{x/hf + C 2 ] / 2 -) } , 

with Cik and C2 being the i nteg ration constants. 
A general solution of Eqs.dl^),(|T^) can be written as 



»=4 

ft V 1 



A-ikfik- 



(20) 



The odd fundamental mode / 4 grows exponentially at 
|x| 3> Iq, and so does the particular solution ft. This 
growth cannot be canceled in Eq.(|20|) by the terms of 
the odd particular solution f^ k . Therefore, one should 
put A^k — 0. The coefficients Ci k , A 2 k and A% k should 
be chosen such that the terms of even particular solu- 
tions f^ k , ft^ proportional to cosh 2 (x/7o) and growing 
exponentially at large x, cancel each other. This gives 

C\k = —fi>A 3 k/ekA2k- 

At distances from the kink, \x\ S> lo, the fundamental 
modes = \x\/x, ft = \x\/Iq — 1, and one can put 
f% = 0. The ratio f k /f k = ^fc/ e fc) an d one can only 
deal with the function f£. Then, in the spatial region 
where l < \x\ < from Eqs.(0)-(|2^) we obtain 



f k + = A lk + e j-A 2k (l - 2^j , x > 
ft = -Ai k + jA 2k (l + , x < 



(21) 



For finding the reflection coefficient we assume that at 
x — > +00 the excitation wavefunction is a superposition 
of incident and reflected waves: 

ft = exp (ikx) + Gexp (— ikx), 

and at x — > — 00 there is only a transmitted wave 

ft = D exp (ikx). 

The reflection coefficient i? = |G| 2 , and the quantity 
|Z?| 2 = 1 — J? is the transmission coefficient. The asymp- 
totic wavefunctions ft (x — > ±00) retain their form in the 
range of distances, Zq <C |x| <C , where they become 



f£ = l + i(l - G)fcr; x > 0, 
/+ = £>(1-Hfcc); x > 0, 



(22) 



The wavefunction (|22|) should coincide with f k (|21|), 
which immediately gives the reflection coefficient 



R = e 2 /V 



(23) 



The reflection coefficient increases with e k and becomes 
of order unity at e k ~ /1. At larger energies i?(fc) strongly 
decreases. Note that for single particles the kink-wise 
density dip is absolutely transparent (R = 0, see fl3|). 



VI. TRANSPORT COEFFICIENTS AND 
LIFE-TIME OF A KINK STATE 

We can now calculate the transport coefficients and 
estimate the life-time of a kink state. At low temper- 
atures T <C /i the main contribution to the transport 
coefficients comes from the scattering of excitations with 
energies e k ~ T. The calculation from Eqs.(^),(||), with 
R(k) (§D, gives 



m !2C(4) r4 

J% TTfl 3 



B 



m 12C(5) T5 



h 



(24) 



The diffusion time tjj (|T^) and the life-time ( |i"3| ) depend 
on the mass M of the kink and, hence, on the transverse 
size of the condensate. For a harmonic transverse con- 
finement with frequency lo±_ the transverse cross section 
S = h/muj±, and from Eqs.(0),(H|) and (§3) we obtain 



T D 1 =24C(4V ± (^ a 3 ) 1 / 2 
t- 1 =2t d x log" 1 [Mfi/mT] 



(25) 
(26) 



where the Riemann (- function £(4) « 1.08. 

In the opposite limiting case, where the trans- 

port coefficients are determined by the scattering of ex- 
citations with energies e k ~ (i. For these excitations the 
reflection coefficient R ~ 1, and we find a ~ mT/Ti and 
£> ~ m/j,T/h. This gives 
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r^ 1 ~ ^(noa 3 ) 1 / 2 -; (27) 
A* 

t- 1 = 2t Z3 1 log" 1 (M/m). (28) 

In the presence of harmonic confining potential in the 
axial direction the Hamiltonian (|5|) of a subsonic kink 
acquires an extra term —Moj 2 q 2 /A, where u> is the ax- 
ial trap frequency ||. This means that in the absence 
of dissipation the kink undergoes oscillations with fre- 
quency Lu/y/2, studied in ^,0]. Importantly, at finite 
temperatures the discrete structure of the spectrum of 
thermal excitations does not manifest itself in their scat- 
tering from the kink, since the level spacing for quasi- 
classical longitudinal (axial) excitations is very close to 
huj/^/2. Hence, the transport coefficients and the diffu- 
sion time td will be the same as in the absence of axial 
confinement. In the regime of friction- induced accelera- 
tion (t > td) the zero-temperature equation of motion 
for the kink Iq] contains an extra term —q/rn' 

q-± + ^-q = 0. (29) 
td 2 

The character of evolution of the kink in a trap de- 
pends on the parameter lotu- For lotd 3> 1 the solu- 
tion of Eq. (|29] ) represents harmonic oscillations with fre- 
quency ui/y/2 and slowly increasing amplitude: q(t) = 
q(0) exp(t/2r£)) sin(wt/v2 + <j)). Hence, the kink reaches 
the border of the condensate and disappears on a time 
scale r = 2td log(Z c /g(0)), where q(0) is the initial am- 
plitude of the oscillations, and l c — (2/i/rraw 2 ) 1 / 2 the ax- 
ial size of the condensate. The amplitude q(0) can be 
found from the condition that in the end of the diffu- 
sion stage of the evolution the kinetic energy of the kink 
Mq 2 /2 ~ Mcj 2 q(0) 2 /4. This gives q(0) ~ {t^/Muj 2 ) 1 / 2 . 
Accordingly, the life-time r of the kink state practically 
coincides with that in a homogeneous condensate and 
follows from Eqs.(|26|),(p8|). This means that near the 
border of the condensate the kink velocity is comparable 
with the velocity of sound. 

For LJT£> <C 1, from Eq.(p9[) we obtain q(t) = 
q(0) exp (t/rp), and the kink acquires the sound velocity 
and disappears before reaching the border of the conden- 
sate. Thus, the dissipative evolution will be the same as 
in the case of an axially homogeneous condensate. 

For Na condensates at densities uq w 3 x 10 14 cm~ 3 
the chemical potential /1 w 200 nK. Then, according to 
Eqs.(fj5l),(|26|), for the radial confinement with tu± w 2kHz 
(satisfying the criterion of dynamical stability for the ax- 
ially Thomas-Fermi condensate, fi < 2.6hu± |t]]) one 
should have T < 50 nK in order to reach the life-time 
of the kink state r > Is. For Rb condensate at the same 
no and wj_ the life-time r > Is requires temperatures 
T < 15 nk. One can think of achieving these conditions 
in an optically confined cigar-shaped condensate, similar 
to the one in the MIT experiment H] . 



VII. CONCLUSIONS 

In conclusion, we have developed a theory of dissipative 
dynamics of a kink state in a finite-temperature Bose- 
condensed gas and found that the dynamics is fundamen- 
tally different from that of vortices H . As the kink has a 
negative mass, due to friction-induced energy losses it ac- 
celerates towards the velocity of sound, and the kink state 
continuously transforms to the ground-state condensate. 
This makes the kink dynamics fast and insensitive to the 
longitudinal (axial) size of the condensate. The fast dissi- 
pative dynamics of kinks is important for understanding 
the picture of relaxation under a rapid quench of strongly 
elongated condensates, where one expects the formation 
of a number of kinks (see, e.g. in Jig]). Importantly, by 
decreasing temperature well below the chemical potential 
one can make r sufficiently large for studying analogies 
between the standing waves of matter and light. 
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